Conversion of ^°K-^^Rb mixtures into stable molecules 
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We study the conversion of ^"^K and *^Rb atoms into stable molecules through the stimulated 
Raman adiabatic passage (STIRAP) in photoassociation assisted with Feshbach resonance. Starting 
with the mean-field Langrange density, we show that the atom-to-molecule conversion efficiency by 
STIRAP aided by Feshbach resonance is much larger than that by bare Feshbach resonance. We 
also study the influence of the population imbalance on the atom-to-molecule conversion. 
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I. INTRODUCTION 

The study on cold atoms is a remarkable research area 
which has been extended from monoatomic systems to 
diatomic systems in recent years [H, 0, H, 0, S S S 111 ■ In 
these experiments, not only degenerate Fermi-Fermi [l], 
0, 0| but also Fermi-Bose 0, |j, H, H, Hi mixtures are 
studied. Principally, the degenerate atoms can be con- 
verted into molecules through the resonant photoassoci- 
ation or magnetoassociation (Fechbash resonance) . Note 
that such a compound molecule created by the Fechbash 
resonance is in a quasibound state and hence energeti- 
cally unstable. Although the quasibound molecules are 
energetically unstable, most of fermionic atoms can be ef- 
ficiently converted into molecules with a long lifetime in 
experiments [HI- However, as far as we know, the 

bare Feshbach resonance can not convert either bosonic 
atoms or Fermi-Bose mixtures into molecules with high 
efficiency. 

To investigate various novel features of cold molecules, 
one must create ground-state molecules with high atom- 
to-molecule conversion efficiency. The bare stimu- 
lated Raman adiabatic passage (STIRAP) [l^, [3 in 
photoassociation is proposed to enhance the atom-to- 
molecule conversion efficiency through avoiding quasi- 
bound molecules' radiation decay. The success of the 
STIRAP technique relies on the existence of the coher- 
ent population trapping (CPT) state i.e., the sys- 
tem can evolve adiabatically in coherent superposition 
of stable states. It is easily satisfied for non-interaction 
systems. However, the inter-particle collisions make it 
difficult for interaction systems to evolve adiabatically, 
which limits a practical efficiency of STIRAP. Compar- 
ing with the bare STIRAP in photoassociation, the stim- 
ulated Raman adiabatic passage aided by magnetoasso- 
ciation is found to be a more efficient technique whose 
conversion efficiency is not limited by the collisions be- 
tween atoms [l^, [l^ [13] • This technique can be applied 
to a monoatomic system or a diatomic system. The free 
atomic states together with the quasibound and ground 
molecular states constitute a four-level system for which 
STIRAP aided by magnetoassociation is applicable. 

In this paper, we consider systems consisting of 
fermionic and bosonic atoms and their compounded 



fermionic molecules through the technique of STIRAP 
aided by magnetoassociation. We take the '"'K-^^Rb mix- 
ture as an example to show that this technique can con- 
vert the two species of atoms into molecules with high 
efficiency. In comparison with STIRAP aided by magne- 
toassociation, we also show that the bare Feshbach reso- 
nance can hardly convert ""^K-^^Rb atoms into molecules. 
The influence of the population imbalance between two 
species on the atom-to-molecule conversion efficiency is 
also studied. In the next section, we present our model. 
In Sec. mil we derive the mean-field dynamical equa- 
tions through the Euler-Lagrange equation then obtain 
the CPT solutions and the corresponding "two-photon" 
resonance condition. In Sec. IIVI we solve the dynami- 
cal equations numerically and discuss the corresponding 
results. Our main conclusions are summarized in Sec. [V] 



II. THEORETICAL MODEL 

We consider a mixture of fermionic and bosonic atoms 
which are coupled to a quasibound molecular state via 
Feshbach resonance. Meanwhile, a laser field drives tran- 
sitions between the quasibound and ground molecular 
states. For convenience, let | /) and | h) stand for the 
ground states of fermionic and bosonic atoms in the open 
channel, and | m) and | g) for the quasibound and ground 
molecular states in the close channel, respectively. The 
state I m) is coupled with states | /) and | h) through 
a magnetic field with coupling strength a' and detuning 
£' . Additionally, the states | m) and | g) are coupled 
with each other through a laser field with the coupling 
strength Vl' and detuning A'. Then the Hamiltonian de- 
scribing the above system in the interaction picture is 
written as. 



H = 



+ ^[*+vI'^^, + H.C.] + (A' + £')*^*3 



(1) 
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where ^'j^) and ^'^j) {h j = f:b,m and g) are the 
annihilation and creation field operators. They obey 
the commutation (+) or anticommutation (-) relations 
i), ^^(r, t)]± = Sij6{r — r') for bosons or fermions, 
respectively. The coefficient Ti is the kinetic energy due 
to the particles' motions, and X[j — X'j^ — 2Trh'^a^j/m^j 
the interaction strength between particles with Oj^- being 
the s-wave scattering length and m^j = m^rrij / [m^ + m^) 
being the reduced mass. The trapping potential term is 
not included in the Hamiltonian as we merely consider a 
uniform system. 

The system we considered includes both bosonic 
and fermionic components. The self-interaction of the 
bosonic component is distinctly different from that of the 
fermionic one. The kinetic energy dominates the intra- 
species interaction for the fermionic component as there 
is no s-wave scattering for two fermions in the same inter- 
nal state. This is in marked contrast to the bosonic com- 
ponent for which the interaction energy dominates the 
kinetic one under most experimental conditions. Thus 
a very good first approximation is neglecting the intra- 
species interaction for fermions but neglecting the kinetic 
energy for bosons. Based on the above consideration, the 
energy density corresponding to the Hamiltonian Eq. ([T]), 
in the Hartree approximation, is given by 

+ ^W^^fi^b + H.C.] ^[rm^, + H.C.] 



'^Kb\^b\ 
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where ipi represents the complex probability amplitude 
of the ith component and A'^ — ?i^(67r^')^/^/2mi. The 
effective self-interaction term Syl^lT/iil^^/^/S related to 
fermions is called Pauli blocking term. One can find that 
the effective self-interaction is in different power of |-0| for 
bosonic and fermionic components. This will induce dis- 
tinct difference between the dynamical equation for the 
bosonic component and that for the fermionic one. 



(see Ref. [l9'|). In the mean-field approximation, the evo- 
lution of each kind of bosons is determined by a single 
equation, however, the evolution of Nf fermions is deter- 
mined by Nf equations. Obviously, it is easy to handle 
the system with bosons but difficult to handle the system 
containing a large number of fermions. Another approach 
is based on the Euler-Lagrange equation with the help of 
the mean-field Lagrangian density [20, 21]. In this case, 
the approximation attributes to a single evolution equa- 
tion for each kind of particles. The difference between the 
mean-field dynamical equations for bosons and that for 
fermions arises from the effective self-interaction term. 

The method adopted in Refs. [13, [ll] is not applicable 
to our system due to the existence of fermions. We study 
the system in terms of mean-field Lagrange density, 



dt 



E. 



(3) 



Substituting the above mean-field Lagrangian density 

1 1 X . dC ^ f dC 

mto the Euler-Lagrange equation —— - d^, 

dip* \d{duipi, 
0, one can get a set of equations for the complex proba- 
bility amplitudes, V'/j i^b^ V'm' and -0^. These equations 
are shown to guarantee the following identities. 



which means that the total numbers of species / and 



b are conserved, i.e., IV'/P 



rif and 



iV'hP + IV'mP + \i>gV = ^b '^it^ ^^'^ ^b being con- 



stants. In the following discussion, we assume there is 
no molecules in the system at the initial time, hence 
and also denote the initial atom densities of the cor- 
responding species. To simplify the calculation, we let 
(pi = '^il \J'nj + n^, then the aforementioned dynamical 
equations for ipi become 

^ Y.^MHs^Af{<p)<Pffi^f + '^^l<p^, 



dt 



/ . '^bi\H'i\ 



a 



III. MEAN-FIELD DYNAMICS AND CPT 
STATE 

The mean-field approach is an effective method to 
solve many-body problems (particularly valid for systems 
with a large number of particles) although the high or- 
der quantum correlations are ignored in this approxima- 
tion 1^ . There are different approaches to get the mean- 
field dynamical equations from Hamiltonian ([1]). A rig- 
orous approach is expanding bosonic and fermionic field 
operators in terms of conventional creation and annihila- 
tion operators, subsequently, substituting them into the 
Heisenberg equations of motion for the field operators 



dt 



E 

n 

~2 



Kni\4'i\'^4>7n + A 



A 



with the conservation relations \(j)f 



ijf'Pb 



(4) 



(1 + S)/2 and 



(1 - (5)/2. Here 



we introduce S = (n^ — n^^)l [n.j to characterize the 
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population imbalance between species / and 5, together 
with further simplified notions: Xij = Kji''^/ + ^b)l^^ 

£ = £ /h, and A = A' /h. The population imbalance 
between two-species atoms is an important parameter 
affecting the feature of the system. For example, the 
population imbalance between fermionic atoms for dif- 
ferent spin state can induce the superfluid to normal 
state phase transition. 

Unlike the monoatomic system Tt'], the densities of 
the fermionic and bosonic atoms may be different, which 
will affect the conventional two-photon resonance con- 
dition. Additionally, the nonlinear terms in our system 
will be changed because of the Pauli exclusive principle 
for fermions. We will see that the STIRAP technique 
aided by Feshbach resonance can convert the two species 
of atoms into molecules with high efficiency, strictly in 
contrast to the situation with Feshbach resonance only. 
We assume that Eqs. (jH) support a CPT steady state with 
(^m = (the validity of this assumption is supported by 
our result obtained in the following). We search steady 



state solutions of Eqs. ([4]) with the help of the following 
trial wave functions. 



0™,, = |0„,,|e^''™.«e-*W+^^.)*, (5) 

where /i^ and /ij, are undetermined parameters. Substi- 
tuting Eqs. ([5]) into Eqs. ^ and taking \(f>m\ = 0, one 
can find a set of solutions, 



,0,2 _ ja^S - O^) + ^ {a^S - + 2a^n^l + S) 



2a2 

= |^0p_5,|^0|2^i + ^ ,.0,2 

together with a restriction condition 



/I ^ - Of ^ "h^ 



(6) 



J 



A = + (A 



(A 



bb 



+ ^fb ^ \g)\ 



r 



(A 



A, 



(7) 



This is a result valid for converting mixtures of fermionic 
and bossonic atoms into fermionic molecules. The pa- 
rameters Xij and Ai can be obtained directly for concrete 
systems. 



IV. NUMERICAL RESULTS FOR CONCRETE 
SYSTEMS 

Now we are in the position to consider a concrete 
system consisting of two-species atoms (saying '^''K and 
*^Rb) and their compounded fermionic molecules. In or- 
der to carry out the numerical calculations, it is nec- 
essary to fix magnitudes of the parameters that affect 
the atom-to-molecule conversion efficiency. As we know, 
£' ~ (2/ie + /^Ar)(^ ~ where /ig and /ijy refer to 
the electron and nucleus magnetic moments, B the mag- 
netic field (the resonances occur at S = Bm) [11] ■ To 
evaluate the parameter a', let us recall the form of the 
s-wave scattering length between species / and b: o-fb = 
a(l — Apcs/{B — B„i)) . Here a is the scattering length far 

off resonance and Apcs — <^'^i^fb^ {'2'!^fi^a{2fie+fJ-pf)) the 
resonance width with m^;, = m^mj,/ (m^ -t-m^) being the 
reduced mass of species / and b. The scattering length a 
for the ■*°K-^^Rb mixture determined through Feshbach 
spectroscopy is about — ISSag with oq being the Bohr ra- 
dius p^, and the expectation value for the width Apcs 
is about -3G \2^. Thus we obtain a' « 9.07 x IQ-^Sj. 



To optimize the conversion efficiency in experiment, 
one can change the detuning £' and the interaction be- 
tween species / and b by varying the magnetic field B. 
Alternatively, one can also fix the magnetic field but vary 
n' to improve the conversion efficiency. In our numer- 
ical calculation, the magnetic field B is fixed to sim- 
plify the experimental procedure. For B = 546. 4G, the 
corresponding s-wave scattering length between species 
/ and b is about 2950ao. If the total atomic density 
+ n^^ is about lO^'^m^'^ and the s-wave scattering 
length between *^Rb atoms is about lOOao, it is easy 
to obtain A^^ = 0.23a, A^^ = 0.0056a, Af = 0.3a, 



An 



A„ 



^fb - 0-23a, ^^bb 
= 0.09a and £ 



-4.1a. Since there are 



no good estimations on molecular scattering lengths, we 
take the interaction strengths involving molecules to be 
zero. The full set of Eqs. ^ are solved numerically by 
choo sing a time-dependent Rabi frequency adopted in 
Ref. [11 



1 — tanh 



t - tr 



(8) 



and A given by Eq. ([7]). Here r2max7 and r are con- 
stants to be determined by the laser coupling of the two 
molecular states. We take fimax — 200a, to = 120/a, 
T = 40/a and assume that there exist no molecules 
in the system at the initial time, i.e., 4'm,g = at 
t = 0. Such an assumption can be realized in exper- 
iment through the following procedure. At the initial 
time, fix the magnetic field at the value far off resonance 



4 



and prepare °K and *'Rb atoms in their ground states 
I /) =1 F = 9/2,mp = -9/2) and | b) =\ 1, 1), respec- 
tively. Turn on the laser pulse and fix its amplitude at 
^^max, and then let the magnetic field be suddenly close 
to the Feshbash resonance point. Some numerical results 
on atom-to-molecule conversion with B = 546. 4G are 
plotted in Fig. [TJ Figure [TJa) exhibits that at 

the initial time, which implies the populations fulfil the 
CPT solutions while a small deviation begins to appear 
at about t = 300 /a. In spite of the small deviation, the 
final conversion efficiency 2\<j)g{t — oo)p is about 60%. 
Note that the Feshbach resonance can not convert ^'^K- 
^^Rb atoms into molecules with such a high efficiency. 
We plot the conversion efficiency as a function of S in 
Fig. [Ijb). Clearly, there is a significant influence of S on 
the conversion efficiency. 
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FIG. 2: (color online) The atom-to-molecule conversion effi- 
ciencies 2\ipg{t = oo)p versus the magnetic field B for differ- 
ent population imbalance S. The parameters except A^j, and 
£ are the same as that in Fig. [l] (a) . 
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FIG. 1: (color online) Panel (a) shows the time dependence 
of the population of particles for S — 0. Panel (b) is the 



conversion efficiency x = 2| 
the population imbalance S. 
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i>g{t = oo)| as the function of 
The parameters are A^j = 0.23, 
X^a =0, Af ^ 0.3, 



bm 



A 



A,„ ^ Ag = 0.09 and £ = -4.1 



bg 



Time is in unit of 1/a and 
all other coefficients are in units of a. 



As we know, except Ayj, and £, the other parameters do 
not change with respect to the magnetic field. Thus the 
magnetic field affects the conversion efficiency through 
the detuning £ and the interaction between species / 
and b. To show the influence of the magnetic field B, 
we plot the atom-to-molecule conversion efficiency ver- 
sus the magnetic field in Fig. [21 One can see that the 
strength of magnetic field can affect the conversion ef- 
ficiency distinctly and the curves are not continuous at 
B = 546. 6G. Such a discontinuity is brought in by the 
divergence of the interaction A^f, at i? = 546. 6G. 

In comparison to the STIRAP aided by Feshbach res- 
onance technique, the results for the bare Feshbach reso- 
nance are plotted in Fig. [3] The results are obtained by 
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FIG. 3: (color online) The time dependence of molecular 
population for the Feshbach resonance. The parameters are 



Aji, = 0.23, 



0.0056, A 



fm 



0, A 



f 



0.3, 



5 = 0, _ 

A,n — 0.09 and £ — —4.1. Time is in unit of 1/a and all other 
coefficients are in units of a. 
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solving the following set of equations: 



V. SUMMARY 



' dt 
' dt 

.d4>m 



J2 XM^bfcbf + Afi^fcbj) 

i—{b,m} 
i={f,b,m} 



rn ) ^ 4'm 



(9) 



Unlike the STIRAP technique, the magnetic field is swept 
according to _B = i?ini + (here Bjni denotes the initial 
value of the magnetic field) in the Feshbach resonance ap- 
proach. The interaction A/b and the detuning £ therefore 
change with time in the calculation procedure if 7 ^ 0. 
From Fig. [3l we find that the Feshbach resonance tech- 
nique can hardly convert ^°K-*^Rb atoms into molecules 
even if the magnetic field is swept slowly. The magnetic 
field is brought from the initial value 551G to the final 
value 547. 5G in 2ms in panel (a), whereas the magnetic 
field is fixed in panel (b) due to 7 = 0. 



We studied the atom-to-molecule conversion efficiency 
for Fermi-Bose mixtures. We took the system consisting 
of ''°K-*^Rb atoms and their compounded molecules as 
an example to have shown that the STIRAP aided by 
Feshbach resonance is an efhcient scheme for one to con- 
vert Fermi-Bose atoms into fermionic molecules. Such a 
scheme is easy to handle in experiments because it is 
no more necessary to sweep the magnetic field across 
the Feshbach resonance. In contrast to the STIRAP 
aided by Feshbach resonance technique, the bare Fer- 
shbach resonance can hardly convert Fermi-Bose atoms 
into molecules. We showed that the population imbal- 
ance 5 between two-species atoms is a conserved quan- 
tity which is an important parameter affecting the fea- 
ture of system. We discussed the influence of 5 on the 
atom-to-molecule conversion in the system that consists 
of Fermi-Bose atoms and their compounded molecules. 
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